arXiv:astro-ph/0512194vl 7 Dec 2005 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN style file v2.2) 


Constraints on Alternatives to Supermassive Black Holes 


M. Coleman Miller^ 

E-mail: miller@astro.umd.edu 

^ University of Maryland, Department of Astronomy, College Park, MD 20742-2421, USA; 
and Goddard Space Flight Center, Greenbelt, MD, USA 


Submitted 2005 October 


ABSTRACT 

Observations of the centers of galaxies continue to evolve, and it is useful to take a 
fresh look at the constraints that exist on alternatives to supermassive black holes 
at their centers. We discuss constraints complementary to those of Maoz (1998) and 
demonstrate that an extremely wide range of other possibilities can be excluded. In 
particular, we present the new argument that for the velocity dispersions inferred 
for many galactic nuclei, even binaries made of point masses cannot stave off core 
collapse because hard binaries are so tight that they merge via emission of gravitational 
radiation before they can engage in three-body or four-body interactions. We also 
show that under these conditions core collapse leads inevitably to runaway growth of 
a central black hole with a significant fraction of the initial mass, regardless of the 
masses of the individual stars. For clusters of noninteracting low-mass objects (from 
low-mass stars to elementary particles), relaxation of stars and compact objects that 
pass inside the dark region will be accelerated by interactions with the dark mass. 
If the dark region is instead a self-supported object such as a fermion ball, then if 
stellar-mass black holes exist they will collide with the object, settle, and consume it. 
The net result is that the keyhole through which alternatives to supermassive black 
holes must pass is substantially smaller and more contrived than it was even a few 
years ago. 
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1 INTRODUCTION 

High-resolution observations of the nuclei of many galaxies 
have revealed large dark masses in small regions. These are 
most naturally interpreted as supermassive black holes, but 
as emphasized by Maoz (1998) it is important to take stock 
of how rigorously we can rule out other possibilities. 

Here we present arguments showing that under ex¬ 
tremely general conditions almost all other options are ruled 
out, further emphasizing that supermassive black holes are 
by far the least exotic and most reasonable explanations for 
the data in many specific sources. In § 2 we lay out our as¬ 
sumptions, making them as conservative as possible so that 
our conclusions are robust. In § 3 we show that for many ob¬ 
served galactic nuclei, binaries are unable to heat the stellar 
distribution effectively because if they are hard then they 
merge quickly via gravitational radiation. This important 
constraint, which depends only on dynamics and not the 
detailed properties of the specific objects, was not presented 
by Maoz (1998) or elsewhere as far as we are aware. In § 4 
we explore the consequences of core collapse and demon¬ 


strate that a very significant mass will inevitably coalesce 
even for point masses. In § 5 we investigate for the first time 
the consequences if stellar-mass black holes exist outside the 
nucleus. We show that enough of them will find their way 
to the center that they will have serious effects on the nu¬ 
clear region, likely consuming a significant amount of mass 
and leading to a supermassive black hole. We discuss the 
consequences of this analysis in § 6. 


2 ASSUMPTIONS AND DYNAMICS 

In the spirit of Maoz (1998), we make a series of conservative 
assumptions to rule out alternatives to supermassive black 
holes. Let us suppose that observations have revealed that 
a mass M is confined within a spherically symmetric region 
whose radius is at most R. We also assume that this mass 
is composed of identical point masses m; the point mass as¬ 
sumption minimises the interaction between the masses, and 
making them identical increases as much as possible the re- 
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laxation time, on which the masses concentrate in the center 
of the distribution and hence increase interaction rates. The 
local two-body relaxation time for a mass m in a region of 
mass density p and velocity dispersion a is 

1 

3 In A G^mp 


where the orbital e nergy per objec t is equal to the kinetic en¬ 
ergy of field stars (IOuinlarJll996a^ . Suppose that the stellar 
velocity dispersion is Vies at the resolution radius Tres for a 
particular galactic nucleus. Then at the hard/soft boundary 
the semimajor axis a is given by 

2Gmla « . (4) 


where InA ~ 10 is the Coulomb logarithm. In general this 
time depends on radius, but note that if p ~ ^.-3/2 
velocity dispersion is dominated by a single large mass, the 
relaxation time is constant with radius. 

Any streaming motion (e.g., rotation) reduces the rel¬ 
ative S£eed_o_and_Jiencereduces the relaxation time (see, 
e.g., iKim. Lee, fc SDurzernll2004l for a numerical treatment 
of a rotating stellar system). Therefore, completely random 
motion leads to the largest timescales. 

For N identical masses in a region whose crossing 
time is fcross, the global r elaxation time is approximately 
iSinnev fc Tremainjll98'^ 

4^ 0.14Ar 4 

trlx ~ ln(0.4JV) ( 2 ] 

« 10® yrMg/®(l MQ/m)(R/l pc)®/® ^ 

where M = lO^MsM©. Both expressions for the relaxation 
time show that for hxed mass density, lower-mass objects 
take longer to alter their distribution, as is expected be¬ 
cause two-body relaxation occurs due to graininess of the 
gravitational potential, which is less when there are more 
objects. 

The more concentrated the initial density distribution 
is, the shorter will be the central rela xation time (see the 
extensive discussion in IOuinlarJll996lJl . To be conservative, 
we therefore assume a relatively flat distribution such as a 
Plummer sphere, in which p oc (1 -I- r®/r®)“®/®, where 
is the core radius. Even for such a distribution, identical 
point masses will undergo core collaps e within a time (see 
discussion in iBinnev fc Tremainelll987^ 

tec ~ 16Aix,h (3) 

where trix.h is the relaxation time at the half-mass radius. 
Note that this is a factor ~ 20 times less than the time 
neede d for the cluster to evaporate iBinnev fc Tremainel 
Core collapse of single objects will formally lead to 
infinite density at the center. In globular clusters and sim¬ 
ilar systems, this is avoided by the intervention of bina¬ 
ries: three-body and four-body scattering can transfer en¬ 
ergy from binaries to the stellar velocity dispersion, heat- 
ing the cluster and stabilising the densi t y at the center (see 
i Gao et ~ ] Il99li : iFregeau et al. ] 120031 : iGiersz fc Spurzeml 
I 2 OOOI for cluster simulations involving primordial binaries). 
As we now show, however, when the velocity dispersion is 
high enough (as it is in many observed galactic nuclei), bi¬ 
naries cannot prevent core collapse. 


3 THE INSUFFICIENCY OF BINARIES 

As shown hrst bv lHeggi3 (^7^, binary-single interactions 
tend to harden hard binaries, and soften soft binaries. Only 
hardening will inject energy into the cluster and slow core 
collapse, hence we only need to consider hard binaries. For 
equal-mass objects the hard/soft boundary is approximately 


Any binary emits gravitational radiation as it orbits. If the 
time for the binary to merge by gravitational wave emission 
is less than the time for the binary to interact with field 
stars, then the binary does not heat the cluster. For a fixed 
semimajor axis, the merger time is maximised for a circular 
orbit, so we assume e = 0 to be conservative. For compar¬ 
ison, if e « 0.7 (the mean for a thermal distribution), the 
merger time is decreased by a factor ~ 10 for fixed a. The 
rate of change in the semimajor axis from gravitational ra¬ 
diation, andmorre^ondlng merger time for a circular orbit, 
is then JPeterslll96^ 

da/dt = —^G®pm®;n/(c®a®) 

Tinerge = a/\da/dt\ = j|g C® (G®771®) (5) 

= |(c/Ures)®(Gm/'!;®es) 

where muin = rni + m 2 is the total mass of the binary, 
p — mim2/7T),bin is the reduced mass, in the second line we 
assume mi = m 2 = m, and in the third line we substitute 
a = 2Gm/i;®es- 

The timescale for a three-body interaction is r 3 _bod = 
l/(?iEu), where n is the number density, v ~ \/2nres is the 
relative speed, and S = irr® [l + 2G(mbin + 77i)/(rpU?es)] is 
the interaction cross section, where Vp is the distance of clos¬ 
est approach. For Vp ~ a and three equal masses, a binary 
at the hard/soft boundary has E ~ dira® « 167rG®m®/u(fes- 
Substituting n = pjm we find 

"Ts-bod « n®es/ (l6v^7rpG®m) . (6) 

The ratio between the merger and three-body timescales is 
then 

I’merge/l'S —bod ~ 44(c/Ures) G pm !■ d 

This ratio needs to exceed unity for the typical binary to 
interact before it merges. Using the average density p « 
p = M/(47ri?®/3) and assuming a roughly constant velocity 
dispersion Vcea = GM/R, we End after some manipulation 
that Tmerge/rs-bod > 1 implies 

m > |(UrBs/c)®/®M 

^2OM0i;f/®3M8 

where Ures = 10®7;res,3 km s“®. A cluster made of any point 
masses lighter than this cannot support itself by binary heat¬ 
ing. 

A loophole might appear to be that when there is bulk 
rotation (and hence a reduced velocity dispersion) or a den¬ 
sity profile in which the relative speed at the center is much 
less than (GM/i?)^/®, binaries wide enough not to merge 
quickly could still heat the distribution. However, suppose 
that a binary has tightened by interactions to the point that 
a = 2Gm/7;®ea, as considered above. Its specific binding en¬ 
ergy is then Gp/(2a) = Gm/Ca) = GM/{8R), because 
w®es = GM/R and thus a = 2R{m/M). Even if the clus¬ 
ter is 100% binaries, the total binding energy liberated by 
hardening is therefore GM®/(8i7). The minimum binding 
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Figure 1. Inferred black hole masses a nd stellar speeds at res- 
olution radius, derived from Table II o f jFerr arese Fordll2Q05^ 
with updates for M31 iBender et al. II2OQ5I) . and the Milky Way 
iGhez et al. ]|2^, which is far to the right of the diagram at 
Vres = 1-2 X 10^ km s“^. The curved lines are labeled by the min¬ 
imum mass of identical point masses such that if they make up 
the dark mass, binaries can in principle heat the system and delay 
core collapse. Several galaxies have Mmin > 100 Mq (the Galaxy 
has Mjnin ^ 400Mq) and hence no reasonable stellar component 
could heat the system. 


energy of a cluster with mass M and outer radius R is ob- 
taiued when all the mass is in a thin spherical shell at radius 
R, in which case the binding energy is GM^/{4R). Even in 
this case, therefore, the maximum effect of binares (prior 
to their reaching the previously considered semimajor axis 
a = 2Gm/v^gg) is to increase the cluster binding energy, and 
hence the cluster radius, by 50%. A smaller binary fraction, a 
more concentrated cluster, or nonzero eccentricities for the 
binaries will all reduce this number. Therefore, if binaries 
that are hard relative to Vres merge quickly by gravitational 
radiation, no possible configuration of velocities or densities 
can allow binaries to stall collapse signihcantly. 

Figure Q plots the black hole mass versus the stellar 
velocity at the resolution radius, along with the minimum 
mass of point masses that would allow binary heating. Sev¬ 
eral galactic nuclei cannot be heated by masses lower than 
100 Mq, including the Galaxy, M87, M31, and NGC 4258. 
If such masses were assembled, the number of objects would 
therefore be small for a given dark mass, which would re¬ 
duce the relaxation time dramatically (see equation 0 and 
would mean that the evaporation time fevap ~ 300 frix would 
be much less than a Hubble time. Therefore, even with an 
implausible collection of > 100 Mq objects in binaries, the 
cluster would still disintegrate rapidly. 


4 CORE COLLAPSE 

If core collapse happens, what is the result? ICohnI (Il980ll 
found that the density profile approaches n oc . For 


ease of calculation, and to be conservative, we will assume a 
shallower prohle n oc r“^, which is appropriate for a singular 
isothermal sphere. In such a prohle, the total mass interior 
to radius r is proportional to r, and the velocity dispersion 
is constant with radius. 

In the high density central regions, even point masses 
can merge because the y emit gravitational radiation. 
lOuinlan fc Shanirol ^1989^ showed that for a relative speed 
V at inhnity between two masses with reduced mass /i and 
total mass mtot, there will be a mutual capture if the peri- 
center distance of approach Vp satishes 

rp < ■ (9) 


For equal masses and v « \/2ures, the cross section for merg¬ 
ing in the gravitationally focused limit is 


27rr p,max(Cmtot / u 



.( 10 ) 


Over a time T, the probability of merger of an average point 
mass is then P = TnSurea. The average number density is 
h = (M/m)/(A tvR^/3). At this density, we Hud after some 
algebra that the probability is 


P 


AT 


M 



V 


3 

res 


GM ■ 


( 11 ) 


With the rough approximation that n « n{rlR)~^ and M(< 
r) « {r/R)M, this implies that the enclosed mass Mmeige 
inside of which the masses merge in time T = lO^Tg Gyr is 




« 5 X 10® MQT^/^{m/l ^ ^ 


or just M if P > 1. The net result is that even for low-mass 
objects, core collapse will lead to the formation of a large 
single mass at the center of the distribution. However, as is 
clear from Equation 0 , if the component masses are small 
enough then the relaxation time is so large that core collapse 
will not occur. We now address this situation. 


5 DYNAMICAL FRICTION AND 
STELLAR-MASS BLACK HOLES 

Suppose that the particles comprising the matter are very 
low-mass indeed, such as elementary particles. Suppose also 
that, like hypothesised dark matter, the particles interact 
neither with themselves nor with ordinary baryonic matter 
in any way but gravitationally. If in some improbable cir¬ 
cumstance the particles have collected in a cluster of total 
mass M and radius R, what will affect them? 

Because the particles have low mass, any more massive 
objects that enter their region will sink to the center via 
dynamical friction. The characteristic time for a mass m to 
sink is (see iBinnev fc Tremainelll98^ for a discussion) 

tdf ~ v%[/ \A-K^\nKG'^fym\ (13) 

where vm is the speed of the massive object. In A is a 
Coulomb logarithm, and 

^ = erf(A) - (14) 

VTT 
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with X = vm/(\/2(t). If vm ~ I'res, then ^ « 0.2. Adopting 
as before Wres ~ {GM/R^^ we find 

tdf ~ 0.2(M/m)(GM/nres) 

« 8 X 10® yrM|(l 

This implies that for systems such as the central region 
of M31, where v^ea ~ 2000 km s“^ and M > 10® Mq 
iBender et al. ll200!^^ . even ordinary stars will sink to the 
center of the mass distribution within a few Gyr, or much 
less if the dark matter is more concentrated. Therefore, all 
ordinary stellar processes that would proceed around a su- 
permassive black hole will also proceed around a concen¬ 
trated region of noninteracting particles, except t hat stars 
inside the region will sink to the center rapidly Isee lOuinlarJ 
Il996hil . Thus if the dynamical friction time at the average 
density p = M/{4'kR^/3) is less than a few Gyr, stars and 
compact objects that enter the region will collide, merge, 
and have prime conditions for forming a large single mass. 

The rate of interactions of stars with the central con¬ 
centrated region is less for smaller regions. Suppose that 
the non-stellar matter is very concentrated, say with a ra¬ 
dius just a few times the radius of a black hole with the 
same mass. Then, the arguments used to estimate rates 
of extreme mass ratio inspirals also apply here. These ar¬ 
guments suggest that stellar-mass black holes will spiral 
into super massive black holes a t a rate not less than ~ 
10~® yr~^ llHils fc Bended llQQ.'j: 1^ gmdsson& 

M^ldarTlscud^^^Gou^]200^]Reitae n20(Bn200^ ~ llvan^ 

200^ ~ lHoDmar^r^Gexandeill2(ffl^ . Therefore, regardless of 
how compactly the dark matter is distributed, if stellar-mass 
black holes exist they will enter the mass distribution in 
much less than a Hubble time. 

The mass accreted by a black hole during inspiral 
is comparatively small. For example, consider a constant- 
density region p = p = M/ with nonrelativistic 

particles moving at an average speed Ures = {GM/R)^^^ 
relative to the black hole. The cross section for absorption 
by a black hole of mass m is T, — {AGm/c?){2Gm/v^g^), so 
during a time tdf the black hole will accrete a mass 


Am = pEUresTDF 

« 0.4 (Ures/c)® m . 


(16) 


This is therefore only a small fraction of the original mass. 
Similarly, if after inspiral the black hole is fixed at the center 
of the mass distribution, it accretes little mass. 

This conclusion changes if the black hole wanders freely 
around the dark matter distribution. This could happen if, 
for example, multiple massive objects enter the dark matter 
region and scatter each other frequently. In this case, for the 
same assumptions as before, the mass accretion rate m = 
pEures becomes 


rh 


Gc2 


(17) 


implying a growth time 


Tg^owth = k{M/m) {GMCjo^) 

« 2 X 10^‘‘ . 


(18) 


This is not constraining on most supermassive black hole 
candidates, but for the Galaxy (M « 4 x 10® M©; 

iGhe^^jgOO^ and Ures ~ 1.2 x 10"^ km s“^ within 45 AU 
iGhe^^i* J20^, the growth time is Tgrowth ~ 4 x 10® yr. 


Radi o observations iReid fc BrunthaleJ 120041 : IShen et al. I 
l2005^ suggest that at least 4 x 10® M© is contained within 
~ 0.5 AU of the position of Sgr A*, which lowers the ac¬ 
cretion time to at most ~ 100 yr. Note that a doubling of 
mass decreases the time to the next doubling by a factor of 
two, so substantial growth results in a runaway. Note also 
that because by assumption the only matter entering the 
black hole does so with prompt infall and without release of 
radiation, the growth is not limited by the Eddington rate. 
If the particles are baryonic or otherwise have a reasonable 
strength of self-interaction, then the accretion rate is greatly 
enhanced, up to a possible Eddington-like maximum. 

Finally, suppose that the non-luminous matter is in 
fact in the form of a star supported by pressure gradients 
rather than by simple motion as we have assumed up to 
this point. An example w ould be fermion balls (see, e.g., 
iTsiklauri fc Viollierlll993l . which are collections of massive 
neutrinos supported by degeneracy pressure. In that case, 
clearly a stellar-mass black hole captured by the star will 
consume matter at the star’s center and remove pressure 
support, leading to rapid destruction of the star. 


6 CONCLUSIONS 

In the past decade, thanks to many observational develop¬ 
ments, the case for supermassive black holes in the centers 
of many galaxies has gone from strong to essentially in¬ 
escapable. We have shown that for many specihc galactic 
nuclei, the observational constraints are strong enough to 
rule out binary heating, hence the relevant evolution time is 
the time to core collapse. This is a factor of ~ 20 less than 
the time to evaporation, which has previously been used as 
the conservative standard for stellar cluster persistence. For 
many individual galactic nuclei, therefore, the combination 
of time to core collapse and lack of binary heating rules out 
dense stellar clusters as an alternate explanation for the in¬ 
ferred dark mass. Specifically, the Galaxy, NGG 4208, and 
M31 have core collapse times <2 Gyr for 0.5 M© objects 
and cannot be stabilised by binaries less than 100 M©. M32 
also has a core collapse time <2 Gyr, but could in principle 
be stabilised by stellar-mass binaries. All other sources cur¬ 
rently have core collapse times >200 Gyr for 0.5 M© objects. 

The only remaining possibilities are concentrated re¬ 
gions of noninteracting low-mass partic les or self-supported 
exoti c objects such as a fermion balls llTsiklauri fc Viollierl 
Il99fjl . Even in this case, we have shown that dynamical 
evolution of the stars and black holes near the centers of 
galaxies will cause multiple stellar-mass black holes to fall 
to the center of the potential, if black holes exist at all. Such 
black holes would consume any high-mass exotic pressure- 
supported objects, and would also accrete a noninteracting 
cluster of particles if allowed to move around freely. There¬ 
fore, the existence of stellar-mass black holes would lead to 
the production of supermassive black holes in many specific 
sources even if the supermassive holes did not form in other 
ways. When combined with the high re dshifts inferred from 
Fe K g lines in some Seyfert galaxies llRevnolds fc NowakI 
|20 o 3), dramatic deviations from standard physics are re¬ 
quired to explain observations in ways not involving black 
holes. 
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